data=xlsread("C:\Users\臻\Documents\WXWork\1688857255593315\Cache\File\2025-04\附件1：123家有信贷记录企业的相关数据.xlsx",'进项发票信息');
result1=zeros(123,11);%记录表
for i=1:123
    %原始数据
    m=find(data(:,1)==1);%企业数据位置
    re=data(m(1,1):max(m),:);%统一公司数据
    re=sortrows(re,[3,4,5]);%按照时间顺序排序
    void=find(re(:,10)==0);%作废发票记录
    void=size(void,1);%作废记录
    valid=size(m,1)-void;%有效记录
    %筛除完作废发票和有效发票中小于0的数据
    mm=find((re(:,10)==0)&(re(:,9)>0));%企业数据位置
    ree=re(mm,:);%筛除完作废发票和有效发票中小于0的数据
    ree=sortrows(ree,[3,4,5]);%按照时间顺序排序
    %计算结果
smonth(1,:)=ree(1,3:4);%发票起始月份
smonth(2,:)=ree(size(ree,1),3:4);%发票终止月份
length=(((smonth(2,1)-2017)*12+smonth(2,2))-((smonth(1,1)-2017)*12+smonth(1,2)))+1;%计算运营时长
up=data(m(1,1):max(m),6);
up=numel(unique(up));
resultl(i,1)=i;
resultl(i,2)=sum(ree(:,7),1);
resultl(i,3)=sum(ree(:,8),1);
resultl(i,4)=sum(ree(:,9),1);
resultl(i,5)=numel(m);
resultl(i,6)=valid;
resultl(i,7)=void;
resultl(i,8)=valid/numel(m);
resultl(i,9)=up;
resultl(i,10)=length;
stdj=std(ree(:,7))./mean(ree(:,7));
resultl(i,11)=stdj;
end
xlswrite("C:\Users\臻\Documents\WXWork\1688857255593315\Cache\File\2025-04\附件1：123家有信贷记录企业的相关数据.xlsx",resultl,"进项结果","A1:k123");
